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Photons, due to the virtually vanishing photon-photon interaction, constitute to very good ap- 
proximation an ideal Bose gas, but owing to the vanishing chemical potential a (free) photon gas does 
not show Bose-Einstein condensation. However, this is not necessarily true for a lower-dimensional 
photon gas. By means of a fluorescence induced thermalization process in an optical microcavity 
one can achieve a thermal photon gas with freely adjustable chemical potential. Experimentally, 
we have observed thermalization and subsequently Bose-Einstein condensation of the photon gas 
at room temperature. In this paper, we give a detailed description of the experiment, which is 
based on a dye-filled optical microcavity, acting as a white-wall box for photons. Thermalization is 
achieved in a photon number-conserving way by photon scattering off the dye molecules, and the 
cavity mirrors both provide an effective photon mass and a confining potential - key prerequisites for 
the Bose-Einstein condensation of photons. The experimental results are in good agreement with 
both a statistical and a simple rate equation model, describing the properties of the thermalized 
photon gas. 
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I. INTRODUCTION 

Photons have played a vital role for the experimental 
realization of Bose-Einstein condensation (BEC) fT in ul- 
tracold atomic gases [5]-[7] . To cool down matter to suffi- 
ciently low temperatures, fundamentally new cooling and 
trapping techniques based on the interaction with light 
had to be developed. These experimental techniques were 
rendered possible only by the invention of the laser and 
even between the birth of the laser and the experimental 
realization of Bose-Einstein condensation in atomic gases 
more than thirty years elapsed. 

Considering their bosonic nature it seems natural to 
ask if photons can undergo a BEC in their own right. 
It is known that a similar phenomenon occurs, when a 
mode in a laser resonator is macroscopically occupied 
above the lasing threshold. The oscillation build-up of 
a laser mode is however not a thermodynamic phase 
transition in a conventional sense. Neither the state of 
the light field is in thermal equilibrium nor the lasing 
threshold can be determined by thermodynamic consid- 
erations. To some extent a laser is even a prime example 
for a system in non-equilibrium: Typically laser gain in 
a medium is only achieved if the transition at the las- 
ing wavelength is inverted, which can only be achieved 
under non-equilibrium conditions using external pump- 
ing [51-[TU|. The question, whether a BEC of photons in 
a strict sense is possible, typically is negated with refer- 
ence to black body radiation |lll IT^ - presumably the 
most omnipresent Bose gas at all [T3|, [T4] . For black 
body radiation not only the spectral distribution of the 
photon energies depends on temperature, but also their 
total number. Temperature and photon number cannot 
be tuned independently. If the temperature of the photon 
gas is lowered then also the photon number will decrease 
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according to the Stefan-Boltzmann law. preventing the 
onset of Bose-Einstein condensation. 



The last years have seen increasing effort in investigat- 
ing equilibrium processes that lead to macroscopically 
occupied photonic states. It is clear that such processes 
have to conserve the photon number. Presumably the 
first process to be considered was Compton scattering 
of x-ray photons in plasmas |15j . Later theoretical work 
considered the use of cavities for photon condensation, 
but did not give a route for reaching the thermaliza- 
tion [16 . In an interesting work, photon-photon scat- 
tering in a nonlinear medium was investigated |17fl26j . 
Here photons are expected to thermalize due to four- 
wave mixing which can be interpreted as collisions be- 
tween two photons - quite similar to binary collisions of 
atoms. However, this ansatz could not be verified ex- 
perimentally. The photon-photon interactions in the ex- 
perimental systems investigated up to now seem to be 
too weak to achieve a thermalization [12]. In other work, 
condensation of classical waves in nonlinear or disordered 
media has been investigated, though the systems are non- 
thermodynamic |27H30| . Experiments on exciton polari- 
ton condensation |31f[55] can also be put into the broader 
context of light condensation. These systems are char- 
acterized by a strong coupling of photonic and exitonic 
degrees of freedom which leads to new energy eigenstates, 
the exciton polaritons. Due to their excitonic part, which 
is coherently coupled to the photonic part, the polaritons 
perform binary collisions, which can lead to a thermal- 
ization of the polariton gas. 



Recently we have observed a photon number conserv- 
ing thermalization process of a two-dimensional photon 
gas in a dye microresonator [3T. Here a thermal con- 
tact between photons and a dye solution working as a 
heat bath is achieved. This leads to a thermalization 
of the transversal photon state. Signatures of this ther- 
malization process are Bose-Einstein distributed photon 
energies and a spatial relaxation of the photon gas in- 
duced by an effective trapping potential. Once the phase 
space density approaches a value of order unity Bose- 
Einstein condensation to the transversal ground mode 
can be observed |3S] . It is one of the purposes of this pa- 
per to give a more detailed account on the experiments 
we have performed so far. In section [H] we describe the 
basic ideas behind the experiment, together with a rate 
equation model describing the thermalization. In section 
|III| the experimental apparatus is described. The experi- 
mental results obtained concerning the thermalization of 
the photon gas are presented in section |TV] followed by 
a detailed account on the experiments on Bose-Einstein 
condensation of photons in section |V] We conclude the 
paper with an outlook on possible future directions for 
the photonic BEG. 



II. EXPERIMENTAL SCHEME 

In our experiment we use a bispherical microresonator 
filled with a drop of a dye solution, see Fig. [l] By 
multiple fluorescence and reabsorption the photon gas 
gets in thermal contact with the dye solution, and takes 
over its temperature, which in our case is room tem- 
perature This thermalization process of the photon gas 
stems from a thermalization process of the rovibronic 
dye state. Due to frequent collisions, on a femtosec- 
ond timescale at room temperature, with the molecules 
of the solvent the rovibronic dye state is permanently 
altered. These processes are fast enough that both ab- 
sorption and emission of photons always take place from 
highly equilibrated rovibronic dye states. As a conse- 
quence the spectral distributions of absorption a{uj) and 
fluorescence f{uj) are linked by the Boltzmann factor 
f[uj)la(uj) oc uj'^ exp {—huj /k-QT). This relation is known 
as the Kennard-Stepanov relation |i39ri42j . Here a{uj) is 
the absorption coefficient of the dye solution, and /(w) is 
the average spectral energy density of an emitted photon 
in free space. By multiple absorption and emission pro- 
cesses this Boltzmann factor is transferred to the spectral 
distribution of the photon gas. This can be seen by inves- 
tigating the corresponding rate equations for absorption 
and emission [13]. More precisely the photon gas ther- 
malizes at the spectral temperature defined by 

where dimensionless profiles, a(w) — Q;(cj)/a(wo) and 
/(w) = f{u})oj~"^/f{(jJo)ojQ^ for an arbitrary ujq, have 
been used [3H]- However, this spectral temperature co- 
incides for a variety of dye molecules with the thermo- 
dynamic temperature of the dye, Tspec(w) ~ T, to good 
approximation. This is especially true for the dyes used 
in the experiment, rhodamine 6G and perylene diimide. 
Regarding the photon number, it is clear that the fiuores- 
cence induced thermalization process conserves the pho- 
ton number on average. A dye molecule can only emit 
an optical photon from its electronically excited state if 
it has absorbed a photon before. A purely thermal ex- 
citation is largely suppressed at room temperature by 
a factor exp[— /iWc/ZcbT] ~ 10~^^ with Huic — 2.1 eV. In 
contrast to a black body radiator, the number of (optical) 
photons in the cavity is not determined by the temper- 
ature of the dye solution. The chemical potential of the 
photon gas thus in general is non-zero. In the experi- 
ment, photons are injected into the cavity by means of 
an external pumping laser, whose power thus determines 
the chemical potential. 

The mirror separation is in the micrometer regime, 
namely Dq ~ 1.46 jim on the optical axis, and is of the 
order of the optical wavelength. Consequently, the free 
spectral range of the resonator, of order ~ 100 nm, be- 
comes comparable to the spectral width of the dye fiuo- 
rescence. The presence of the mirrors modifies the spon- 
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Figure 1: Experimental scheme, (a) Resonant modes of the 
microcavity (top) and relative absorption a{u) as well as flu- 
orescence strength /(i/) for rhodamine 6G (bottom). The 
height of the bars indicates the degeneracy of the photon 
energy. Inside the microcavity, to a good approximation 
only fluorescence photons with a longitudinal wave number 
of g = 7 (black bars) exist, which causes this degree of free- 
dom to be frozen out making the photon gas effectively two- 
dimensional, (b) Schematic set-up of the dye-flUed microcav- 
ity. The curvature of the mirrors provides a trapping poten- 
tial in the transversal plane, as indicated on the left. Photons 
with high transversal excitation numbers (green) are emitted 
at a larger angle to the optical axis, whereas transversally less 
excited radiation (yellow) exhibits a lower spatial divergence. 



taneous emission from the excited molecules inside the 
cavity such that predominantly photons of a certain lon- 
gitudinal mode number q are found to be emitted, typi- 
cally g = 7 in our case. By this, one degree of freedom is 
frozen out and the photon gas becomes two-dimensional, 
as only the two transversal mode numbers are variable. 
Moreover, this establishes a ground state for the pho- 
tons with non- vanishing energy - the transversal ground 
mode TEMgoo- In thermal equilibrium, the transversal 
excitations are expected to be Bose-Einstein distributed. 
Consequently, for increasing temperature of the dye so- 
lution, the average angle of propagation of the photons 
with respect to the optical axis will increase, while for 
low temperatures the photons are expected to propagate 
nearly parallel to the optical axis. 



A. Cavity photon dispersion 

Formally, the photon gas can be mapped onto a two- 
dimensional, ideal Bose gas in a harmonic trapping po- 
tential. This can be seen by investigating the energy- 
momentum relation. The energy of a photon as a func- 
tion of longitudinal wave number and transversal 

wavenumber kr = 



fc^ -|- fc2 jg given by 

Eph = —y/k^ + k^, 
n 



(2) 



where n is the index of refraction of the medium. Due 
to the mirror curvature the boundary conditions for the 
photon modes (in z-direction) depend on the distance 
from the optical axis r = |r|. For simplicity these bound- 
ary conditions are taken to be metallic. For the longitu- 
diiral wavenumber k^ = kz{r) one can make the ansatz 



where 



k^(r) =^ qii/D{r), 



D{r) =Do- 2{R - \J B? - r'^) 



(3) 



(4) 



is the mirror separatioir at a distance r from the optical 
axis and R is the radius of curvature of the mirrors. To 
incorporate a possible nonlinear photon interaction we 
moreover allow the index of refraction of the medium to 
vary in space n = n{f) = uq + An^, with An^ = n2/(r), 
in the presence of the light field. Here uq and n2 are the 
linear and nonlinear indices of refraction of the medium, 
respectively, and 1(f) is the intensity. A non-zero value 
for 712 can arise e.g. from a Kerr-nonlinearity or from 
thermo-optical effects. With that we derive 
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where in the last step both the paraxial approximation 
(r <^ R, kr -€1 kz) and Art^ ^ no has been applied. For 
an effective photon mass defined by 
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equation ([s]) can be written as 
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When the interaction term is omitted, i.e. for n2 = 0, this 
is the energy-momentum relation of a non-relativistic. 
massive particle in two dimensions, bound by a harmonic 
trapping potential. The standard quantization procedure 
then leads to the well known energy spectrum of the bi- 
spherical resonator |44j, with the energy levels equidis- 
tantly separated by Ml and the degeneracy increasing 
linearly. After rescaling the energy such that the lowest 
harmonic oscillator mode has zero energy, and defining 
the transversal photon energy u = iJph — rnpi^c^ /uq — Ml 
the degeneracy can be written as 



g{u) = 2{u/hVl + 1). 



(9) 



The prefactor 2 here goes back to the two-fold polariza- 
tion degeneracy. 



B. Statistical theory of the photon gas 

The thermodynamics of such a system is well known: 
a harmonically confined two-dimensional Bose gas shows 
a phase transition to a BEC at low temperatures or high 
particle numbers [45 -48_. The average occupation num- 
ber of a certain transversal energy u is given by the Bose- 
Einstein distribution 
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(10) 



where the chemical potential of the photons is implicitly 
defined by 



(11) 



with A^ph as the total (average) photon number. With 
that, the number of photons in the transversally excited 
resonator states is seen to saturate at a critical particle 
number given by 



Nr. 



kBT 



(12) 



If the total photon number is increased beyond this 
bound, a macroscopic fraction is expected to condense 
into the transversal ground state. In our experiments, 
the trapping frequency is of order Q, ~ 10^^ Hz, which 
is roughly 9 orders of magnitude higher than for exper- 
iments with ultracold atomic gases. At room tempera- 
ture, T — 300 K, this gives a critical photon number of 
magnitude Nc — 10^, which is experimentally feasible. In 
principle, it would also be possible to lower the temper- 
ature of the dye solution below the critical temperature. 



^/?,Ml 
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Trfce^-o y Do R 



ph 



(13) 



to trigger the condensation for a fixed photon number. 
However, this is experimentally much more demanding. 



Strictly speaking, phase transitions can only occur in the 
thermodynamic limit. From eq. 13 we see that this ther- 
modynamic limit has to be Aph — ^ oo, i? — > oo with 
Np\i/B. = const, as only this retains T^. A limit invoking 
Dq — > oo is unphysical because the two-dimensionality of 
the photon gas will be lost. Other than in the confined 
case considered here, there is no thermodynamic limit 
for a macroscopic ground state occupation at finite tem- 
peratures for the homogeneous two-dimensional Bose gas 
|45fl48] . Correspondingly, such a system is always multi- 
mode in the thermodynamic limit. 

In contrast to polariton experiments |31fl36| , no strong 
atom-light coupling takes place in our experiment. The 
frequent collisions with the molecules of the solvent leads 
to a rapid decoherence of the dye state, i.e. the decoher- 
ence rate is bigger than the Rabi frequency. This prevents 
a coherent energy exchange between photonic and molec- 
ular excitations |49l [50] . Consequently, the eigenstates of 
the system are best described by the uncoupled photonic 
and molecular degrees of freedom. At any given point 
in time a certain number of photons, A'ph, and excited 
dye molecules, Aexci populate the cavity. In a station- 
ary state, i.e. no net fiux, the ratio of these numbers will 
be linked by the ratio of the lifetimes of both species, 
Acxc/A^ph — Texc/Tph- The lifetime of a photon Tph here 
denotes the average time between emission and reabsorp- 
tion in the medium. For typical experimental parameters 
Tph is of magnitude 10 ps. The lifetime of the excited dye 
molecules t^xc is well known to get modified by the pres- 
ence of the cavity with respect to its free space value 
[STHMl- Its exact value has a complicated dependence 
on the detailed experimental conditions. Moreover, it 
will be infiuenced by the presence of other photons due 
to stimulated emission. Nevertheless, we expect it to 
be of order 1 ns for the used optical powers, and thus 
it remains clearly bigger than the photon lifetime un- 
der typical experimental conditions [Ml 155]. As a result, 
the number of excited molecules will also be much bigger 
than the number of photons. This situation is best de- 
scribed by a grandcanonical particle exchange, where the 
photon gas exchanges particles with the reservoir of ex- 
cited dye molecules. The question regarding the suitable 
statistical ensemble is important because the fiuctuations 
of the photon number - and with that the coherence prop- 
erties of the light - will strongly depend on the statistical 
ensemble, even in the thermodynamical limit 



C. Rate equation model 

Semiclassical laser theory builds upon rate equation 
models, where laser action is predicted when the gain 
from stimulated emission exceeds the total loss [TU]. In 
an analogous way we use a multimode theory based 
on rate equations to describe the thermalization in the 
here investigated dye-filled optical microcavity. The dye 
molecule is modelled by a Jablonski energy level scheme 
[63] commonly used in the Kennard-Stepanov theory of 
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Figure 2: (a) Energy level scheme of a dye molecule with 
an electronic ground state So and an excited state S\ re- 
spectively, where "S" denotes singlet, which are separated 
by an energy fiwo- Photons of a frequency uj are emitted and 
absorbed from within the rovibrational substates E and E' 
respectively, (b) Chemical reaction type exchange between 
different particle species, namely ground state dye molecules, 
photons and excited dye molecules. 



absorption and emission |39fH^ [M] . The present calcu- 
lation neglects cavity specific modifications of the light- 
matter interaction, as well as fluctuations of the particle 
number, both issues that will be incorporated in a sub- 
sequent publication [65 . 

Consider a system with one ground and one excited 
electronic state, each with a manifold of rovibrational 
substates, see Fig. [2^. Let Swq denote the energy splitting 
between the lowest rovibronic states (zero-phonon line). 
Furthermore, E and E' are considered to be rovibrational 
energies of ground and excited state and q{E), g'{E') de- 
note the corresponding rovibronic densities of states. The 
distribution functions p{E), as well as p'{E'), are propor- 
tional to the probability of finding a rovibronic state with 
energy E and E'; for example p{E) — exp{—E/kQT). 
With this, we expect the following total rates for spon- 
taneous emission, stimulated emission and absorption of 
a photon at frequency lu: 



R. 



spoil 



(w) = J dE' ne{E') A {E',uj) , (14a) 



W = JdEg{E)p{E) and W = J dE' g' {E')p'{E'). This 
treatment assumes that cavity mirror losses and the cou- 
pling to unconfined optical modes are small, i.e. the rates 
associated with the corresponding loss processes are small 



compared to the processes accounted for in (14a|-(14c 
Similarly, the rate for pumping, which compensates for 
these losses, will then be small. These conditions will be 
valid if a photon collides several times with dye molecules 
in the cavity before being lost, as can be achieved by re- 
absorption, i.e. by recapturing fluorescence. 

The Einstein coefficients are not independent, their 
connection is given by 

A{E', uj) dE' = hujD{Lj)B'{E', uj) dE' , (16) 

(A-B relation), with D{u!) as the spectral mode density 
per unit volume, and by 

g'{E')B'{E',uj)dE' = g{E) B {E , u) dE , (17) 

which follows from a detailed balance argument |64| . If 
we assume that the rovibrational dye state is in ther- 
mal equilibrium in both the electronic ground and ex- 
cited states manifold respectively, which for many dyes 
is a reasonable assumption given that the typical time for 
a collision with solvent molecules is in the femtosecond 
regime at room temperature, while the electronic life- 
times are typically of the order of a nanosecond [63 , we 
can apply Boltzmann distributions for p{E) and p'(E'). 
From the rate equations (14a|-(fT4c| the net variation 



of the photon number n{uj) at frequency uj is 



dn(uj) 
dt 



= i?spon('^) + i?stim(w) - i?abs(w) (18a) 

j dE g{E)e-'^/^^^B{E,Lo) x 
^A^,e-'"'(^-^°)/'==^^[feji:>(cj) + u{uj)] 



1 

W 



NgU{Uj) 



(18b) 



i?stim(w) = j dE' n,{E')B'{E',u)u{uj) , (14b) 

dE ng{E)B {E , w) u{u) , (14c) ^ and ^ to combine all integrals 



where in the second step we have used the energy con- 
servation relation hu + E — hujQ + E' and the equations 



where u(w) is the spectral energy density per unit volume 
of the radiation fleld. A, B and B' are the corresponding 
Einstein coefficients, and the occupation densities ne{E') 
and ng{E) of the excited and ground state populations 
are given by 



In a steady state, dn(w)/dt — 0, the term in squared 
brackets must vanish. This yields the equilibrium photon 
energy distribution 



n{uj) 



Ji(u}-uja)/k-BT _ ^ ' 



(19) 



ne{E') 
ng{E) 



g'{E')p'{E') 
w 

g{E)p{E) 
w 



(15a) 
(15b) 



Here Ng and N^. denote the number of molecules in 
ground and excited states in the optically active vol- 
ume. Further, W and W' are statistical weights given by 



where additionally n{uj) = J dV and J dV D{uj) = 
g{uj), see eq. (|9|, were employed. In chemical equilibrium, 
the chemical potentials of all species are linked by /iph + 
/ig = fie (Fig- ^p)j where /iph is the chemical potential 
of the photons and /ig, /jg are the chemical potentials 
of excited and ground state dye molecules, respectively. 
From this one can easily show [65 that the fugacity of 
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the photons is determined by the excitation level in the 
dye medium 



NgW 



(20) 



or, if we use a renormalized chemical potential fi — /iph - 
huJr. 



NgW 



(21) 



Thus, in accordance with the usual convention, the pho- 
ton number approaches infinity, Npx^ — > oo, for a vanish- 
ing chemical potential, /i — ?> with T = const. When 
we combine the equations ( 19 1 and (21 ), we arrive at the 



Bose-Einstein distribution, eq. (10) 



As it is well known from laser theory |10| . the rate 
equations ( 14a )-( 14c I do not allow for an inversion of 



a specific transition, and one therefore has to include a 
pumping term to the equations to model lasing. How- 
ever, the rate equations do predict thermally distributed 
photon energies for certain conditions, as we have shown. 
This includes that the statistical weight of the resonator 
ground mode can become macroscopic. Put differently, 
'single mode '-behaviour can arise even in thermal equi- 
librium and without inversion. This point is however 
only revealed in a statistical multi-mode treatment and 
is missed in the usual laser theory. That the photon 
condensate can arise in a non-inverted medium can also 
be understood by the capability of our system to reuse 
spontaneous emission. In a macroscopic laser reabsorp- 
tion nearly inevitably leads to photon loss due to spon- 
taneous emission into non-confined modes. However, if 
spontaneous emission is directed into the cavity, where it 
gets reabsorbed, the excitation will not be lost. There- 
fore spontaneous emission in general is not a loss channel 
in our system. Recapturing of fluorescence (or photon re- 
cycling) is a mechanism also known for "zero threshold" 
microlasers |66 -68 . In those works one tries to capture 
spontaneous emission into a single cavity mode, while 
in the present scheme it is sufficient if the spontaneous 
photons are directed into one of the (many) transversal 
modes. The thermalization process will select the res- 
onator ground state to have the biggest statistical weight, 
which can even be macroscopic, when condensation sets 
in. 



III. EXPERIMENTAL SETUP 

The experimental setup, as shown in Fig.|3] consists of 
a microresonator, a pump source and an analyzing sec- 
tion. The microresonator is formed by two highly reffect- 
ing, spherically curved, dielectric mirrors which are typ- 
ically used for cavity ring-down spectroscopy. Their re- 
ffectivity is bigger than 0.99997 for the wavelength region 
(500 — 590) nm, with a reffectivity maximum of 0.999985 
at A w 535 nm. The latter corresponds to a finesse of 



F w 200000 for the empty cavity. Two different radii of 
curvatures are available, i? = Im and i? = 6m, respec- 
tively. To allow for a mirror separation of a few halfwaves 
on the optical axis despite the curvature, one of the mir- 
rors is cut to a surface area of 1 mmx 1 mm. The edges are 
carefully polished to minimize possible stress from con- 
tact points between the mirrors. To verify that the mir- 
ror's reflectivity remains unaffected by the cutting and 
subsequent polishing procedure cavity ring-down mea- 
surements have been performed. 

The microresonator is pumped at an angle of 45° with 
respect to the optical axis. At this angle and for a certain 
polarisation the mirrors show their first reflectivity min- 
imum allowing a pump beam transmission of over 80%. 
To pump the resonator in this way is advantageous be- 
cause resonance conditions and mode matching do not 
have to be considered. The pump beam can be placed 
at an arbitrary position inside the resonator plane inde- 
pendent of the cavity resonances. One of the mirrors is 
mounted on a translation stage with a piezoelectric actu- 
ator to allow fine tuning of the mirror separation, which 
can be computer-controlled to achieve an electronic sta- 
bilization. For that the fiuorescence spectra emitted from 
the cavity are captured by a spectrometer and are used 
as an input for the stabilization protocol that controls 
the voltage driving the piezoelectric actuator. The band- 
width of this control loop is sufficient to compensate ther- 
mally induced fluctuations of the cavity length. Faster 
fluctuations are reduced by mechanically decoupling the 
cavity from the optical table and by the damping caused 
by the dye solution between the mirrors. 

The used dyes axe either rhodamine 6G or perylene- 
diimid (PDI), whose absorption and fluorescence spectra 
are shown in Fig. |4] The dyes show a quantum efficiency 
close to unity in the used spectral region, $r6G ~ 0.95 
and $pDi ~ 0.97 ^D], and have spectral temper- 
atures close to the thermodynamic temperature of the 
surrounding solvent, rspcc(w) — T. For rhodamine we ei- 
ther use methanol or ethylene glycol as solvent, for PDI 
aceton is used. Repeated filtering of the solutions is nec- 
essary to remove unsolved dye particles and other con- 
taminations. This filtering has proven to be crucial for 
the reliability of the experiment. To achieve reabsorption 
of the microcavity radiation fairly high concentrated dye 
solutions have to be employed. In principle high dye con- 
centrations can lead to radiationless deactivation. How- 
ever, the concentration used here, 1.5 x 10~^Mol/l for 
rhodamine 6G, is still one order of magnitude lower than 
the concentration where fiuorescence quenching notably 
sets in i71n73j . As pump source we use a frequency dou- 
bled Nd:YAG laser at a wavelength of 532 nm which is 
switched by two acousto-optic modulators operating in 
series to ensure a high on/off contrast. The analysis is 
performed by a spectrometer and a camera. For that 
the emitted radiation is split by a non-polarizing beam 
splitter. One part of the beam is mapped by an imaging 
system onto the chip of a color CCD camera, where a 
magnified real image of the photon gas is generated. The 
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Figure 3: Schematic set-up of the microcavity experiment. The microresonator is pumped at an angle of 45° to the optical 
axis. The radiation emerging from the microresonator is examined both spatially and spectrally. 
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Figure 4: Relative strengths of absorption Q(A)/amax and 
fluorescence /(A) //max for the dyes (a) rhodamine 6G and 
(b) perylene diimide. 
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other part of the beam is directed into a spectrometer. 
We use either a commercially available one, with opti- 
cal resolution ~ 2nm, or a self-built spectrometer. Care 
has to be taken when the spectrometer has an entry slit, 
as it is the case for the commercial spectrometer, where 
the width is 15 pm. Due to their bigger mode diameter 
and higher divergence the coupling of higher transversal 
modes through the slit typically will be much less effi- 
cient than, for example, for the ground mode TEMqOQ. 
Consequently, the entry slit will not be neutral in color 
anymore. One way to overcome this problem is to use 
a diffusing plate in the light path that cancels out the 
correlation between position, angle and color. A more 
advantageous way is to omit the entry slit alltogether 
and to directly pass the light onto a diffractive grating. 
This is possible because the emitted radiation from the 
microcavity is only weakly divergent and can be suffi- 
ciently coUimated before being spectrally decomposed. 
This scheme is realized in our self-built spectrometer, us- 
ing a rotating diffractive grating with 2400 lines/mm and 
a photomultiplier tube as detector. However, both meth- 
ods described above deliver largely equivalent spectra. 



Figure 5: Spectral distribution of the radiation emitted 
from the microcavity for different mirror separations, D « 
{37, 20, 11, 4.7} /im, corresponding to different cavity orders 
q ~ {175,95,50,22}. (Rhodamine 6G in methanol, g = 
1.5 X 10"'' Mol/1, mirror curvatures _Ri = Im and _R2 = 6m) 



IV. THERMALIZATION OF THE 
TRANSVERSAL PHOTON STATE 

A. Spectral and spatial photon distribution 

Figure [5] shows typical spectral distributions of the 
emitted microresonator light for decreasing mirror sep- 
aration (from upper left to lower right). By investigating 
the free spectral range the longitudinal mode number of 
the photons can be determined. For very small mirror 
separations a regime is reached where the photon gas 
gets two-dimensional, i.e. where only photons of a sin- 
gle longitudinal mode number are emitted. The shortest 
mirror separation that can be reached in our experiment 
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Figure 6: Transmission of the pump beam (squares) as a 
function of the longitudinal mode number q. The solid line 
corresponds to the Lambert-Beer law I{q)/Io = exp(— a(g — 
go)-^c/2no) with the fitting parameter qq = 4.68 ± 0.17. This 
measurement shows that the effective end nodes of the op- 
tical resonance are located inside the mirrors. (Rhodamine 
6G in methanol, £i 2 x 10"^ Mol/1, mirror curvatures 
iii = -R2 = 1 m) 



corresponds to g = 7 halfwaves. The limiting factor that 
prevents a further reduction of q is the penetration of 
the light field into the mirror material. This can be seen 
by determining the transmission coefficient of the pump 
light, which is illustrated in Fig. [6j For this measure- 
ment very high concentrated dye solutions are used (rho- 
damine 6G at 2 X 10~^Mol/l). As expected, the trans- 
mission of the pump light shows an exponential decay 
with increasing q (Lambert-Beer law). However, the de- 
cay curve extrapolates to full transmission at a value of 
qo = 4.68 ± 0.17. Apparently the thickness -Ddyc of the 
dye film is not directly proportional to q but shows an 
offset, following D^yeiq) oc (q — qo). For example, the dye 
film thickness at g = 7 is just over two halfwaves. The 
effective end nodes of the light field are located about 
2.3 halfwaves (corresponding to 400 — 450 nm) inside the 
mirrors. Longitudinal mode numbers g < 5 are therefore 
not possible with the used mirrors. 

The true ground state of the resonator is the halfwave 
resonance TEMioq. The freezing out of the longitudi- 
nal mode number at mirror separations corresponding 
to ? <^ 7 is not obvious at first sight. Fluorescence 
into TEMgniii modes, or into transversally low excited 
TEMgmn modes are energetically excluded due to the fi- 
nite spectral width of the dye emission. But fiuorescence 
into transversally high excited TEMgmn modes is ener- 
getically possible. However, due to their much higher 
mode volume and consequently smaller overlap with an 
emitting point dipole, decays into higher TEMgmn modes 
are expected to be much slower than emission processes 
into the energetically equivalent but transversally low ex- 
cited TEMymn modes. This is indeed also seen in the 
experimental spectra (Fig. 5|, where a preference for low 
transversal modes is already seen for rather large mirror 
separations. Thus mode hopping to longitudinally lower 
excited modes is not found to be a relevant loss mecha- 
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Figure 7: Spectral distribution of the radiation emitted from 
the microcavity far below the critical photon number at 
T = 300 K und T = 365 K (circles). The spectra are in good 
agreement with the Boltzmann-distributed photon energies 
(lines). For comparison a, T — 300 K Boltzmann-distribution 
is additionally plotted in the bottom graph (dashed line). 
(Rhodamine 6G in ethylene glycol, g — 5x 10~^ Mol/1, mirror 
curvatures Ri — R2 = 1 m, q = 7) 
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Figure 8: Normalized spectral distribution of the radiation 
emitted from the microcavity (circles) for 4 different cut- 
off wavelengths Ac « {570, 590, 615, 625} nm. In addition, 
Boltzmann-distributed photon energies for T = 300 K are 
plotted (lines). (Experimental parameters as in Fig. [7| 



nism for the two-dimensional photon gas. 

Measurements of the spectral photon distribution in 
the microresonator at room temperature, T = 300 K, 
and for a heated resonator setup at T = 365 K can be 
seen in Fig. [7] The power of the transmitted light is 
^out = (50 ± 5) nW. From this value an average photon 
number of iV = 60 ± 10 inside the resonator can be de- 
rived. A numerical solution of eq. ( [TT| then leads to a 
chemical potential of At/fceT = -6.76 ±0.17 (T = 300 K) 
and n/kBT = -7.16 ± 0.17 (T = 365 K). These mea- 
surements are thus performed far from the phase tran- 
sition that is expected to set in for n getting close to 
zero. Consequently the term '-1' in the denominator 



of the Bose-Einstein distribution, eq. (10 1, can be ne 



glected and the spectral distribution of the photons is 
expected to be Boltzmann-Hke. Besides some deviations 
around 532 nm caused by scattered pump light the mea- 
sured spectra show a good agreement with the theoreti- 
cal expectations both for T — 300 K and for T = 365 K, 
respectively. This can be interpreted as evidence for a 
thermalization process of the transversal photon state. 

As stated before a thermalization process is expected 
to occur only if a thermal contact between photon gas 
and dye solution is achieved by reabsorption. To ver- 
ify this experimentally we varied the cut-off wavelength 
Ac = ^-Kcjuj^ by tuning the piezo voltage, as is shown in 
Fig. [8] For smaller Ac the absorption coefficient of the 
dye solution gets bigger and consequently the emitted 
photons are reabsorbed faster. Thus a sufficient thermal 
contact is expected to remain in this case. The upper 
curves in Fig. [s] confirm this (note the logarithmic in- 
tensity scale). The thermal drop-off of the transversal 
excitations is indeed visible. This does not hold anymore 
if Ac is moved to a larger wavelength, see the lower two 
spectra in Fig. |8] Here the absorption coefficient is not 
big enough to achieve reabsorption and consequently the 
spectra differ quite strongly from a Boltzmann distribu- 
tion. A similar observation can also be made by varying 
the dye concentration. For low dye concentrations the 
spectra are only partially thermalized. 

Not only the spectral but also the spatial photon dis- 
tribution shows the expected characteristics of a thermal 
gas. In Fig. [9^ a picture of the photon gas taken by a 
color CCD camera is shown. For that the photon gas 
was mapped with a lens system as a magnified real im- 
age onto the CCD chip. Transversally low excited pho- 
tons (yellow) are seen in the direct vicinity of the trap 
centre while the highly excited photons (green) oscillate 
stronger in the transversal plane. From this picture the 
spatial distribution of the photons, here a cut along an 
axis through the centre, can be determined, see Fig. |9]3. 
Additionally the theoretical thermal expectation value, 
which far below criticality is given by a Gauss distribu- 
tion, is plotted here (solid line). The good agreement 
between the two curves again confirms the thermal char- 
acter of the photon gas. This is however not unexpected 
as the spectral distribution of the photons has already 
revealed Boltzmann-like statistical weights and spectral 
and spatial distributions can be mapped onto each other 
in a thermalized system. 



B. Redistribution of photons 

To further characterize the thermalization process we 
investigated the influence of the pump beam position. 
For that the pump beam is displaced with respect to the 
optical axis. Due to the thermalization process we nev- 
ertheless expect the photons to accumulate in the trap 
centre where their potential energy is minimized. The 
result of such a measurement can be found in Fig. |10[ 




Figure 9: (a) Picture of the photon gas (real image onto the 
sensor of a color CCD camera). Low-energy photons (yel- 
low) are emitted from the direct vicinity of the trap centre, 
while the emission of higher-energetic photons (green) occurs 
outside of the centre, (b) Photon distribution along an axis 
intersecting the trap centre (circles), extracted from Fig.[9|a). 
The theoretical distribution function (line) is based on a ther- 
mal averaging over all resonator modes and it corresponds to a 
Gaussian distribution. (Experimental parameters as in Fig.[7| 
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Figure 10: Measured distance of the intensity maximum of 
the emitted fluorescence from the trap centre, |a;max|, as a 
function of the position of the excitation spot, Xexc- For a 
cut-off wavelength of Ac = 620 nm, corresponding to weak 
reabsorption, the fluorescence exactly follows the pump spot 
(circles connected by dashed line), for Ac = 580 nm, corre- 
sponding to strong reabsorption, the photons concentrate to 
the trap centre, x^^^ = 0, provided the pump spot is displaced 
less than 60 iim from the trap centre (squares connected by a 
solid line). (PDI in acetone, q — 0.75 g/1, mirror curvatures 
i?i = i?2 = 1 m, g = 7) 



In this figure the distance of the fluorescence maximum 
from the optical axis as a function of the pump position 
is shown. The pump spot is stepwise moved along a line 
through the centre and at every position a;exc the dis- 
tance between the point of brightest fluorescence and the 
optical axis l^maxl is measured. The diameter (fwhm) of 
the pump spot is roughly 25]im. These measurements 
are performed for two different cut-off wavelengths Ac. 
At Ac ~ 620 nm reabsorption of the fluorescence light is 
nearly absent due to the much smaller absorption coef- 
ficient in this spectral regime and the decreasing mirror 
reflectivity. Consequently the photons leave the cavity 
in the direct vicinity of the pump spot and the opposite 
turning point, i.e. pump spot and absolute value of the 
position of the intensity maximum are nearly identical. 
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For Ac ~ 580 nm reabsorption is strong enough to trigger 
the thermahzation process. Here the majority of photons 
can accumulate in the potential minimum before they get 
lost. However, this only holds if the position of the pump 
beam is not too far away from the optical axis. This max- 
imum distance is determined to be roughly ~ 60 pm for 
the given experimental conditions. If the pump spot is 
moved beyond that radius, the photon gas does not reach 
its equilibrium spatial distribution anymore. The loss of 
spatial thermahzation is accompanied by a non-thermal 
spectral distribution. We interpret this observed spatial 
redistribution effect as further evidence for the presence 
of a thermahzation process. 

The reason for the incompleteness of the thcrmaliza- 
tion process for stronger displaced pump spots is at- 
tributed to the finite number of emission-reabsorption 
cycles a photon can undergo on average before it is lost. 
For larger displacements this number apparently is in- 
sufficient to reach the equilibrium spatial distribution 
of the photon gas. We expect the most important loss 
channel to be the coupling to unconfined modes. The 
mirrors used in the experiments show a drop-off in re- 
flectivity for angles larger than 45°. At every emis- 
sion process therefore a non-negligible probability ex- 
ists that a photon is scattered out of the cavity. Mir- 
rors with a full three dimensional bandgap would cer- 
tainly improve on this [75] . A second loss chan- 
nel is radiationless deactivation. The average number 
of emission-reabsorption cycles is bounded by the quan- 
tum efficiency smaller than unity. This average is given 
by J2n=o - = - For rhodamine 6G 

('I'r6G = 0.95) it corresponds to 19 fluorescence pro- 
cesses, for PDI ($pDi = 0.97) it is 32. Both numbers 
seem to be rather big and we therefore do not expect the 
flnite quantum efficiency to be the dominant loss mech- 
anism. This is especially true because the quantum effi- 
ciencies could even improve compared to the free space 
values we have stated above. Such an effect is expected 
if the lifetime of the excited molecules is shortened due 
to modified spontaneous emission within the cavity. 



V. BOSE-EINSTEIN CONDENSATION 

A. Spectral and spatial photon distribution 

We now describe the behaviour of the photon gas for 
higher intracavity powers. For these experiments the res- 
onator is pumped with powers of Ppump ^ 100 mW. How- 
ever during the one-time pass of the pump beam through 
the cavity only a small fraction of the light is actually ab- 
sorbed in the dye film - less than 1 mW, the remainder 
is simply transmitted. To prevent excessive population 
of dye triplet states and heat development, the pump 
light is chopped into rectangular pulses of 0.5 jis dura- 
tion. The duty cycle is typically 1:16000, i.e. every pulse 
is followed by a dark phase of 8 ms duration. For pulse 
durations longer than 0.5 jis the intracavity power is ob- 
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Figure 11: Spectral photon distribution for increasing intra- 
cavity power. The optical powers are normalized to the exper- 
imentally determined critical power Pc.exp = (1.55 ± 0.60) W, 
which corresponds to a critical photon number of Nc = (6.3 ± 
2.4) X 10* (Rhodamine 6G in methanol, £> = 1.5 x 10"^ Mol/1, 
mirror curvatures Ri — R2 = Im, q = 7, pulse duration 
0.5 ps). The inset features theoretical spectra based on a 
Bose-Einstein distribution of the transversal excitations. 



served to decay. We attribute this to the population of 
triplet states. This time scale is indeed also known for dye 
lasers as long as no triplet quenchers are used [76 . Since 
the pulses are at least two orders of magnitude longer 
than the lifetime of the excited state dye molecules and 
roughly four orders of magnitude longer than the life- 
time of the photons, the experimental conditions can be 
considered quasi-static. 

The experimentally observed spectral distributions of 
the cavity photons for increasing intracavity power is 
shown in Fig. [TT] In the inset the theoretical spectra 
based on Bose-Einstein distributed transversal excita- 
tions at 300 K are depicted. The circulating light power 
inside the cavity is determined by measuring the trans- 
mitted power through one of the mirrors and using the 
transmission coefficient as the proportionality factor. As 
observed before, the spectral distribution is Boltzmann- 
like for low photon numbers. However for higher photon 
numbers the maximum of the distribution slightly shifts 
to larger wavelengths accompanied by a spectral narrow- 
ing. Beyond a certain critical value, a spectrally sharp 
peak near the cut-off appears. Due to the insufficient res- 
olution of the spectrometer of sa 2nm, it is not possible 
to decide from the spectra which mode is massively pop- 
ulated. Clear-cut evidence for the massive population of 
the ground mode is however provided by the spatial dis- 
tribution shown in Fig. [12] The central bright spot visible 
has a diameter (fwhm) of (14 ± 2) jim. This corresponds 
well to the expected diameter of the transversal ground 
state dxEMqoo = 2-y//i ln2/mphri — 12.2 jim for the given 
experimental parameters. The experimentally deter- 
mined critical light power inside the cavity is Pc,exp ~ 
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Figure 12: Picture of the photon gas (real image onto the sensor of a color CCD camera). Low-energy photons (yellow) 
are emitted from the direct vicinity of the trap centre, while the emission of higher-energetic photons (green) occurs outside 
of the centre, (b) Photon distributions along an axis intersecting the trap centre for increasing order parameter Nq/N « 
{0%, 0%, 1%, 3%, 10%, 25%} (circles), (c) Interference pattern obtained by superimposing two partial beams of the condensate 
using a M ich elson- interferometer. Here the path difference of the two interferometer arms is 15 mm. (experimental parameters 
as in Fig. 11 1 



(1.55 ± 0.60) W. If this value is normalized by the power 
per photon, Pph — /iWo/Tit, where r^t is the round trip 
time, one obtains a critical photon number of iVc_oxp — 
-Pc.expT'rt /^wo ~ (6.3 ±2.4) X 10'*. Within the quoted un- 
certainty, the value is in agreement with the theoretically 
expected number of = {n'^ /3){kBT/hn)^ ~ 77000 for 
the given parameter and thus confirms that the observed 
condensation phenomenon occurs at a particle number 
which is characteristic for a BEC. We have repeated the 
measurement for a different dye (perylene-diimide (PDI) , 
0.75 g/1 in acetone). Since the only property of the dye 
that effects the condensation process should be its tem- 
perature, we do not expect the critical photon number to 
differ from the rhodamine case. This is indeed confirmed 
experimentally; within the given uncertainties no devia- 
tion of Pc,oxp or -^c,exp is observed. There is however a 
difference regarding the pump power that is needed to 
achieve the critical photon number inside the resonator. 
The reason for this is the lower absorption coefficient of 
PDI at the pump wavelength of 532 nm. This is roughly 
a factor of 3 smaller than for rhodamine. To reach the 
same amount of absorbed pump power one therefore has 
to pump harder. 

The relatively large uncertainties regarding Pc.cxp smd 
Nc,exp stem from a combination of mainly systematic er- 
ror sources, especially uncertainties in the calibration of 
the power meter and in the transmission coefficient of 
the mirrors, t = (2.5 ± 0.4) x 10~^. Moreover, since the 
cw-power coupled out at the phase transition is only of 
order ss 1 nW, care has to be taken that the power mea- 
surement is not falsified by scattered pump light or back- 
ground fluorescence from the edges of the mirrors. Here a 
spatial filterering has proved to be helpful. However, one 
has to ensure that especially the higher transversal modes 
are not blocked by this filtering. At last it is also neces- 
sary that the spatial pump distribution, here a Gaussian 
ground mode with beam waist 2wo ~ 100 pm centered 



around the optical axis, differs not too strongly from the 
equilibrium distribution of molecular excitations. As it 
already has been shown in the previous section, the pho- 
ton gas does not reach its equilibrium state if this spatial 
mismatch gets too big. In this case a position dependence 
of the chemical potential remains that is not completely 
removed by the thermalization process before photons 
leave the cavity. Such a non-equilibrium situation will 
also influence the critical photon number. 

To obtain the spatial photon distribution, pictures 
of the photon gas for various pump levels have been 
recorded and analyzed quantitatively. The intensity pro- 
files (cut along an axis through the trap center) shown 
in Fig. |12[ 3 have been normalized to equal area and are 
shifted in the y-direction for the sake of visibility. These 
profiles show an increasing spatial concentration of the 
photons in the trap centre prior to the phase transition 
- a feature, which is in good agreement with the theoret- 
ically calculated intensity profiles. However, there is an 
interesting deviation of the experimental profiles from the 
theoretical ones; the experiments reveal an increase of the 
ground mode diameter with the ground state occupancy. 
This is not expected for an ideal gas of photons and thus 
indicates the presence of a repelling photon-photon inter- 
action induced by the medium. Such an interaction could 
e.g. arise from a Kerr-nonlinearity of the dye molecules. 
However, in our case it seems to stem mainly from a 
thermo-optical effect, further investigations are presented 
in section 5.5. Another indication of the selfinteraction is 
the non-perfect saturation of the thermal modes, which 
can be attributed to a change in the effective trapping 
potential for higher occupation numbers |77j . 

Preliminary experiments regarding the coherence of 
the condensate have been performed. The light beam 
coupled out of the resonator is sent through a Michel- 
son interferometer with a pathway difference of roughly 
15 mm, see Fig. |12| ;. For this measurement the interfer- 
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ometer is slightly misaligned such that the two interfering 
beams do not propagate coUinear. The resulting interfer- 
ence pattern demonstrates that the first order coherence 
of the condensate at least extends into the centimeter 
regime. However, a more thorough characterization of 
the coherence length is subject to ongoing experiments 
and will be presented in a subsequent publication. 



B. Scaling of the critical power with the resonator 
geometry 

The critical circulating power as a function of the ge- 
ometry parameters of the resonator can be casted in the 
form 



(22) 
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One expects a linear increase of with increasing ra- 
dius of curvature R and no dependence on the longi- 
tudinal mode number g - at least as long as the two- 
dimensionality of the photon gas holds. To test the de- 
pendence on R, mirrors with two different curvatures are 
available, i?i = Im and i?2 = 6m, allowing three differ- 
ent mirror combinations. The combination i?ieft = 1 m 
and i?right = 6m is equivalent to a resonator built from 
two identical mirrors with a radius given by the h arn ionic 



mean R = 2 {R^ 



Rr, 



1.71 m. In Fig. 



13 i the 



experimentally determined values of the critical power as 
a function of the effective radius R are shown together 
with the theoretical expectation given by eq. (22 1. The 



data points coincide well with the theoretic prediction 
- both the linear scaling and the absolute values. Seen in 
another way, the procedure of increasing both the radius 
of curvature R and the photon number N but preserving 
the ratio R/N is observed to retain the critical tempera- 
ture of Tc = 300 K. Thus, if this procedure was continued 
indefinitely, one could reach the thermodynamic limit for 
condensation at a finite temperature. In a further mea- 
surement we have investigated the dependence of the crit- 
ical power Pc on the longitudinal mode number q. For 
that, q is increased stepwise without varying the cut-off 
wavelength Ac = 585 nm, see Fig. 13 d. As expected, q 
has no significant influence on Pc within the given exper- 
imental uncertainties. Both measurements again confirm 
that the observed condensation occurs at the BEG criti- 
cality condition. 

Interestingly, transitions to lower longitudinal modes 
were not observed to be significant, which we attribute to 
the vanishing spatial overlap of the modes. The observed 
spectra indicate that the photon gas remains to good ap- 
proximation two-dimensional up to a cavity length of at 
least 16 halfwaves. The pump power that is necessary 
to achieve the critical photon number shows a strong de- 
pendence on the mirror separation, see upper graph in 
Fig. If the longitudinal mode number is doubled 

from g = 8 to g = 16 by increasing the mirror separa- 
tion, only one third of the initial pump power is necessary 
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Figure 13: (a) Critical circulating intracavity power as a func- 
tion of the effective mirror curvature R. (Rhodamine 6G 
in methanol, g — 1.5 x 10"'^ Mol/1, q = 7, pulse duration 
0.5 lis) (b) Pump power at the phase transition (top) and 
critical intracavity power (bottom), respectively, as a func- 
tion of the cavity order q. (Rhodamine 6G in methanol, 
g = 1.5 X 10~^ Mol/1, mirror curvatures _Ri = _R2 = Im, 
pulse duration 0.5 ps) 



to reach criticality. Thus the system shows a behaviour 
contrary to threshold behaviour reported for dye |66} [78] 
or semiconductor microlasers |67| . where an increase of 
the threshold pump power with the cavity length has 
been observed. The reason for the here observed decrease 
of the pump power is an increase of the absorbance for 
thicker dye films. Since the observed critical circulating 
power Pc is independent of the cavity length, the ab- 
sorbed pump power at criticality should also be roughly 
constant. Thus in the weak absorption limit, as in our 
case, one expects to see a reciprocal dependence of the 
pump power Ppump on the dye film thickness Ddye{q), 
as this conserves the amount of absorbed pump power. 
Thus we expect Ppump oc (<? — 5o)~^j with the penetra- 
tion depth into the mirrors qo. The corresponding fit to 
the experimental data is shown in the upper graph of 
Fig. |13| 3. For the parameter qg one obtains a value of 
go = 4.77 ± 0.25, in good agreement with the value of 
qa = 4.68 ± 0.17 derived in section 4.1 by directly mea- 
suring the absorbance. This confirms that the necessary 
absorbed pump power to trigger criticality remains con- 
stant for various mirror separations. The given values of 
the critical pump power can also be converted to pump 
intensities, /pump,c = ^pump,c/7''WQ, where wq ~ 50 jam is 
the pump beam radius. One obtains values roughly in 
the regime (0.6 - 2.2) kW/cm , where the given values of 
P ^ 



pump.c 



in Fig. 13 3 have to be corrected by the mirror 



transmission coefficient, which is of order 50%. These 
values are roughly two orders of magnitude smaller than 
in macroscopic dye lasers, where typically pump intensi- 
ties of order « lOOkW/cm^ are required \29j 80) . 



C. Condensation by spatial relaxation 

As it was already demonstrated in section 4.2, the ther- 
malization process is accompanied by a spatial redistribu- 
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Figure 14: Intensity distribution of the photon gas along an 
axis intersecting the trap centre for varied cut-off wavelengths 
Ac (left edge). The pump beam (dashed line) is located out- 
side the trap centre and its position as well as its power are 
kept fixed during the measurement. The top curve shows a 
photon gas precisely at the phase transition with a ground 
state population Nq/N < 1%. (Rhodamine 6G in methanol, 
g — 1.5 X 10~^ Mol/1, mirror curvatures iii = 7i2 = 1 m, 
q = 7, pulse duration 0.5 jis) 



tion of the photons towards the trap centre. This redistri- 
bution can also be observed for higher photon numbers. 
It can even provide a sufficiently high density at the trap 
centre to trigger the condensation. To demonstrate this, 
we have performed the following experiment: the dye so- 
lution is pumped SOjtm away from the optical axis with 
a pump beam of roughly SSjim (fwhm) diameter. Both 
its position and its power are kept fixed in the course of 
the experiment. By fine tuning the cavity cut-off Ac, we 
can vary the amount of reabsorption. 

The lowest profile in Fig. [14] shows the intensity of the 
photon gas along a line through both trap centre and 
pump position for a cavity cut-off of Ac — 610 nm. Addi- 
tionally the profile of the pump spot alone is displayed, 
where thereto one of the cavity mirrors was removed. For 
this value of Ac the degree of reabsorbtion is low and the 
thermalization process is suppressed. The observed spa- 
tial distribution of the photon gas roughly follows the 
pump beam intensity distribution. If Ac is decreased, 
the amount of reabsorption is increased and the thermal 
contact between photon gas and dye solution is gradually 
established. The intensity profiles get distributed more 
symmetrically around the trap centre. For a cavity cut- 
off of Ac = 570 nm the photon density in the centre even 
gets sufficiently high to trigger criticality. This is both 
indicated by a small bright spot occurring at the posi- 
tion of TEMoo-mode and the corresponding cusp in the 
intensity profile (upper curve in Fig. 
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condensate fraction that is achieved in this way is rather 
low, No/N < 1%. 

An important conclusion from this measurement is 
that the condensation is triggered not directly by the 
pumping. The intensity of the pump beam at the position 
of the TEMoo-niode is almost negligible. We attribute 



this to the multiple emission-reabsorption cycles, the in- 
duced thermalization process provides a sufficient photon 
density for the condensation. Such an effect is unknown 
in lasers but has already been observed in the context of 
polariton condensation[32j . As in section 4.2, a complete 
spatial relaxation of the photon gas is achieved if the 
pump spot is not moved beyond a distance of ~ 50]tm 
from the trap centre. Again this can be attributed to the 
finite number of emission-reabsorption cycles a photon 
can undergo before it is lost. 



D. Reabsorption cycles 

Perhaps the most decisive quantity regarding the de- 
gree of thermalization is the average number of reabsorp- 
tions of a photon before it is lost from the resonator. 
The observed spatial redistribution of the photons to- 
wards the trap centre already tells us qualitatively that 
multiple reabsorptions take place in the experiment. A 
more quantitative measure is experimentally obtainable 
by a time-resolved measurement, for example of the over- 
all decay time of the microresonator fluorescence. Such a 
measurement is under current preparation. Nevertheless, 
already an overall input-output analysis of the data pre- 
sented in Fig. |13| 3 con&ms that the experiments are per- 
formed in a regime of multiple reabsorptions, and gives 
an upper bound for the number of reabsorptions per pho- 
ton. 

From the general conservation of excitations we con- 
clude that the average numbers of excited molecules iVcxc 
and the number of photons Npi^ are linked by 



(23) 



where Tph is the average time between emission and reab- 
sorption of a photon, and Texc the lifetime of the molec- 
ular electronic excitation in the resonator. Further, /3 is 
the probability that an excited molecule emits into a con- 
fined cavity mode rather than decaying radiationless or 
into an unconfined mode. With this equation we implic- 
itly make the approximation that losses predominantly 
occur for conversions from molecular excitations to cav- 
ity photons but not for the inverse direction, which seems 
to be justified in our system. iVoxc is also linked to the 
absorbed pumping power Pp^abs by 



Pr 



p,abs 

fp,abs 
hVr, 



rcxc(l + /? + /3' + ...) 
1 



(24) 



where hv-p is the energy of a pump photon and thus 
Pp.a.hs/hvp is the absorbed pump photon rate. The frac- 
tion of excitations that were directly created by the pump 
beam corresponds to the term '1' of the sum, the fraction 
of excitations that were reabsorbed once corresponds to 
the term '/?', and so on. We here assume that the storage 
time in the resonator is dominated by the lifetime of the 
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molecular excitation, Tcxc ^ ''ph- By the combination of 



eq. ( 23 I und eq. ( 24 1 the probability (3 can be derived 



/3 = 



1 



1 



Tph 



(25) 



and with that, the average number of reabsorptions per 
photon fire = Si^Lo n/3^{l — /?) is written as: 



ph 



(26) 



The lifetime Tph is experimentally controlled by the dye 
concentration and is also strongly dependent on the 
cavity cut-off Ac- For the experimental parameters of 
Fig. 13 i we expect a value of Tph = (21 ± 6) ps. The ab- 
sorbed pump power is estimated to be Pp^abs = (0.65 ± 
0.10) mW. With that we obtain n,.o = 3.8 ± 2.5. The 
error of this estimate is rather big, nevertheless fiio is 
clearly found to be of magnitude 1. The quantitative 
analysis of the particle flux thus confirms that the ex- 
periment is carried out in the regime of multiple scat- 
tering. However, it delivers an upper bound of n^e < 6. 
At first sight this number seems to be relatively small 
compared to the numerous collisions taking place in an 
atomic gas. Regarding the degree of thermalization this 
might seem unfavourable. But one has to consider that, 
in contrast to atomic collisions, nearly all correlations 
between absorbed and emitted photon states are anni- 
hilated already by a single absorption-emission process 
(Kasha's rule[63 | IST ] ) - with the only exception of a nec- 
essary spatial overlap between both photon states. Thus, 
the contact with a heat bath, as it takes place in our pho- 
ton gas experiment, is a much stronger thermalization 
process compared to two-body collisions in an atomic gas. 
Experimentally, we have strong evidence for the thermal- 
ization to operate both properly in the spatial and the 
spectral regime, see chapter [TV| 



E. Selfinteraction of the light condensate 



123) 



The observed spatial intensity distributions (Fig. 
reveal an increase of the condensate diameter wit : 
creasing occupation level. The diameter (fwhm) of the 
transversal ground state as function of N^/Nph is given 
explicitly in Fig. |15| 3. Slightly above criticality the mea- 
sured diameter of (14 ± 2) jim is in good agreement with 
the expected ground state diameter of 



(27) 



for non-interacting photons, yielding do — 14.7 jim for 
the parameters of Fig. [12] For higher ground state occu- 
pation iVo/A^ph this does not hold anymore. This obser- 
vation suggests the presence of a repelling photon-photon 
interaction mediated by the dye solution, which can for- 
mally be described by a non-zero nonlinear index of re- 
fraction, 712 7^ 0, see eq. (|8]). In our case the dominant 



contribution to n2 seems to be given by a thermo-optical 
effect. The massively occupied ground state should lead 
to a small heat gradient in the dye solution due to nonra- 
diative decays. This heat gradient then causes a gradient 
of the refractive index, which decreases the optical dis- 
tance between the mirrors on the optical axis. Effectively, 
this can be interpreted as a local distortion of the mir- 
rors, see Fig. [15^, which is indeed expected to lead to 
a broadening of the ground state mode - in qualitative 
agreement with the experimental observation. 

The index gradient necessary to notably broaden the 
condensate peak can be estimated by comparing he 
transversal ground state energy, HQ, with the interaction 
term E-mt of eq. (|8| 



E\nt = 



(28) 



with An ~ n2/(r). They become equal for a gradi- 
ent of An ~ — fifino/niphC^ which in our experiment 
is roughly An ~ —0.9 x 10^"*. Such a change in the 
refractive index is caused by a temperature gradient of 
AT = {dn/dT)-^An ~ 0.17K, where for methanol we 
used dn/dT ~ -4.86 x lO""^ K"^ [Sg, and thus we expect 
thermal effects to play an important role. In principle, 
higher orders of polarizability as e.g. the optical Kerr- 
effect, contribute to the nonlinear index of refraction. 
This is however expected to be a much weaker process: 
for an index change of magnitude An — 10^*, an inten- 
sity of / ~ An/nj^*"^"^-* ~ IGW/cm^ is necessary, which 
is 4 to 5 orders of magnitude higher than the intensity 
building up in the microresonator. Here a nonlinear index 
of refraction of magnitude nj^''"-' ~ — 1 x 10~^^cm^/W 
was assumed |55] . 

The interaction term, eq. (28 1, can be cast into a more 



familiar form by introducing the condensate wavefunc- 
tion tpo{f) via /(r) = {mphC^ /nlrj-t) No\ipo{r)\'^ , where 
Tit — hqn^/mphC^ is the resonator round trip time. Then 



eq. (28) reads as 



^2 



Tlph 



(29) 



with the dimensionless interaction parameter g given by 



mj^i^c^na 



(30) 



With the interaction term in the form of eq. ( 29 ) the sim- 
ilarity of the photon energy given by eq. ([s]) to the Gross- 
Pitaevskii (GPE) equation becomes even more apparent. 
To determine g in our experiment we measure the occu- 
pation level A^oZ-ZVph at which the condensate diameter is 
doubled compared to the non-interacting case, eq. (27 1. 
This is achieved roughly at Nq/N ~ 0.25, which corre- 
sponds to an occupation number of Nq ~ 40000. From 
a numeric solution of the 2d GPE we obtain a doubling 
of the diameter for an interaction strength of gNo — 30, 
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Figure 15: (a) Schematic illustration of the optical self- 
interaction. The massively occupied ground state is expected 
to locally heat up the dye solution, reducing the index of re- 
fraction and, therefore, shortens the optical distance between 
the mirrors. This is effectively equivalent to a local deforma- 
tion of the mirrors (which is shown exaggerated in the figure), 
resulting in a repulsive interaction, (b) Diameter (fwhm) of 
the fundamental transversal mode (TEMqo) as a function of 
the condensate fraction (circles). The solid line indicates the 
theoretically predicted mode diameter for an interaction pa- 
rameter of ^ = 7.5 X 10"* expected from the Gross-Pitaevskii 
equation, the dashed line gives the diameter expected from a 
Thomas-Fermi approximation, (experimental parameters as 
in Fig. 
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giving the dimensionless interaction paramete r in our ex- 
periment, g ~ 30/40000 = 7(3) x 10"". Fig. [l5]3 shows 
the experimentally determined data points together with 
a curve based on the numerical solution to the 2d GPE 
(solid line) where we used g — 7.5 x 10~^ as an input pa- 
rameter. Apparently the experimental data can be mod- 
elled satisfactorily by this. Additionally, in the figure we 
have included the condensate diameter (fwhm) according 
to a Thomas-Fermi approximation (dashed line), given 

by 



jTF 



2^h/V^mpi,n{gNo)K (31) 



which is supposed to yield approximatively correct results 
for higher condensate fractions, where the transversal ki- 
netic energy of the photons can indeed be neglected. 
Compared to the interaction parameters reported for 



two-dimensional atomic Bose gases, which are in the 
regime g = 10~^ — 10~^ [HI], the interactions in the 
photon gas are at least one order of magnitude smaller. 
Thus the photon gas is closer to an ideal Bose gas, and 
there is no indication that the observed phase transition 
differs from the BEC scenario. This would be expected 
for stronger interacting gases, where the interactions give 
rise to a Kosterlitz-Thouless (KT) type of phase [84 -86]. 
An indicator for this would be the loss of long range spa- 
tial order. In initial experiments using a shearing inter- 
ferometer, see Fig. [12] we have tested for the first-order 
coherence of the condensate by bringing different spatial 
parts of the condensate to overlap. The observed inter- 
ference patterns however did not show any signatures of 
a spatially varying phase. This is in agreement with the- 
oretical expectations, predicting that long range order is 
lost only for higher values of the dimensionless interac- 
tion parameter |84J- 



VI. CONCLUSIONS 

We have investigated thermodynamic properties of a 
two-dimensional photon gas in a dye-filled optical micro- 
cavity. In initial experiments, we have investigated spec- 
tral and spatial properties of the photon gas below criti- 
cality in detail, and obtained evidence for a thermalized 
two-dimensional photon gas with non- vanishing chemical 
potential. When the photon number is increased above 
a critical value, Bose-Einstein condensation of the pho- 
ton gas is observed. Experimental signatures are Bose- 
Einstein distributed transversal excitations including a 
macroscopic population of the transversal ground mode, 
the condensation setting in at the expected critical pho- 
ton number and a spatial relaxation process into the trap 
centre that can lead to condensation even for a displaced 
pump spot. 

The hitherto performed experiments were focussed on 
the thermodynamics of the photon gas, a more detailed 
investigation of the condensate properties is currently 
performed. Initial interferometric measurements of the 
condensate peak have already confirmed the first-order 
coherence of the condensate, with a lower bound for the 
coherence length of the order of centimeters. We also 
plan to investigate the second-order coherence jHT] US] , 
which allows to monitor the unusually large particle fluc- 
tuations expected from a grand-canonical ensemble. This 
would allow to further distinguish the radiation emit- 
ted by the equilibrium photon BEC from nonequilibrium 
laser radiation. 

On a longer timescale, an experiment based on pho- 
tonic crystal mirrors with a three-dimensional bandgap 
|74|[75] should prove advantageous, as losses due to spon- 
taneous emission out of the cavity should be reduced. 
This would allow to reduce the amount of pumping to 
maintain the average photon number in the resonator. 
The full bandgap mirrors together with a condensed mat- 
ter system, e.g. dye molecules in a polymer matrix, could 
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even prove useful for technical applications, as it would 
allow the construction of fluorescence collectors [89, , that 
are more efficient than currently used ones. A tantalizing 
perspective includes the possible realization of coherent 
light sources in the UV-regime. Other than in a laser, no 
inversion is required for the realization of a photon BEG, 
and spontaneous emission is recaptured |66H68| . 

On the fundamental research side, we expect the pho- 
tonic gas within a 'white-wall' box to be a versatile tool. 
For the ideal Bose gas one does not expect superfluidity, 
in contrast to the interacting case. Therefore it would 
be interesting to observe the superfluid properties of the 
condensate as a function of the dimensionless interaction 
parameter g. A possible route for the tuning of g would 
be either the tuning of the photon mass by moving the ex- 
periment to the UV spectral region, the reduction of the 
resonator length by reducing the penetration of the elec- 
tric field into the mirrors, and by increasing the nonlinear 
coefficient rt2. We expect that the interaction parameter 
of currently g Ri 7 x 10^* can be increased by at least 
an order of magnitude, thus allowing to investigate the 
superfluid properties over a wide range of parameters. 



We are aware that the character of the interaction needs 
to be investigated further. In a different line of experi- 
ments, the usage of a second 'heating' dye molecule with 
low quantum efficiency, whose absorption is spectrally 
shifted away from the observed spectral region, should 
allow to thermooptically change the refractive index by 
means of a heating laser. This would allow the modelling 
of e.g. periodic lattices, and also more complex potentials 
for the photon gas, by simply imaging the corresponding 
pattern onto the 'heating' dye molecules. 
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